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Abstract 

In this paper we present the equations of motion of a moving screw dislocation 
in the framework of the translation gauge theory of dislocations. In the gauge field 
theoretical formulation, a dislocation is a massive gauge field. We calculate the 
gauge field theoretical solutions of a uniformly moving screw dislocation. We give 
the subsonic and supersonic solutions. Thus, supersonic dislocations are not forbid- 
den from the field theoretical point of view. We show that the elastic divergences 
at the dislocation core are removed. We also discuss the Mach cones produced by 
supersonic screw dislocations. 
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1 Introduction 

In the dynamics of dis locations it i s usually assumed that the screw disloc ation possesses 



a Lorentz symmetry (iFrankl . 1 19491 ; iHirth fc Lothd . Il968l : iGiintherl . Il988l ) in contrast to 



the edge dislocation, which does not have a Lorentz-type symmetry because of two char- 
acteristic velocities, namely the velocities of transversal and longitudinal waves, entering 
the field equations. By means of such a Lorentz transformation, the screw dislocation 
can be transformed into a steady-state one like in Maxwell's theory of electromagnetic 
fields. The question arises if the 'classical' dynamics of dislocations reproduces the cor- 
rect behavior of dislocations. Usually, the conventional theories of dislocation dynamics 
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have several drawbacks (e.g. inertial effects are missing, singularities). Unlike special 
relativity where the speed of light is an upper limit, in elastodynamics the speed of sound 
is not a limit velocity. Shock waves can move faster than the velocity of sound and 
form Mach cones. Similar predictions that dislocations can move faster than the speed 



e.g.. Eshelbv (1956); Weertman ( 


1967h: Giinthei 


(1968. 


19691): Weertman & Weertman 


(1980): Callias & Markenscofll ( 


1980h: Markenscoff & HuaneJ (2008 


)) and computer sim- 


ulations (Gumbsch & Gaol. 19991: 


Koizumi et al. 


20021: 


Li & Shi. 


2002: Tsuzuki et al. 



2008). Supersonic dislocations have been recently observed in plasma crystals (INosenko et al 



20071 ). 



Another question comes up whether dislocations cause symmetric or asymmetric force 
stresses. In classical dislocation theories the elastic distortion tensor is asymmetric and 
the force stress tensor is assumed to be symmetr ic. Due to this asymmetry between the 
distortion and the stress tensors, IKronerl (119761 ) has introduced a modulus of r o tation 
for t he antisymmet ric part of the force stress tensor. In addition, iHehl fc Kroner! (1l965t ) 
and IKronerl (119681 ) argued th at dislocation s prod uce moment stresses as response and 



antisymmetric force stresses. iPan fc Chen! (11994] ) calculated the antisymmetric s tress 



corres ponding to the body couple using nonlinear and linear continuum mechanics. iGao 
( 119991 ) has developed an asymmetric theory of nonlocal elasticity based on an atomic 
lattice model and the three-dimensional quasicontinuum field theory. He pointed out 
that both strain and local rotation should be regarded as basic variables of geometric 
deformation and he has shown that the local rotation makes a very important contribution 
to the internal energy. For isotropic materials, the antisymmetric stress comes from 
the long range property of interaction of atoms in the metal ma teria l s and non-uniform 
distribution of atomic forces in the more microscopic structures (jGaol . Il999l ). 

A promising and straightforward candidate for an improved dy namical theory of dislo- 
cations is the so-cal l ed translational gauge theory of dislocatio ns (IKadic fc Edelenl . Il983 ; 
Edelen fc Lagoudad . Il988l ; lLazar fc Anastassiadid . 120081 . [2009j). In the gauge theory of 
defects, dislocations arise naturally as a consequence of broken translational symmetry 
and their existence is not required to be postulated a priori. Moreover, such a theory uses 
the field the oretical framework which is well accepted in theoretical physics. 

Recently, ISharma fc Zhangl (120061 ) have applied the gauge theory of dislocations to a 
uniformly moving screw dislocation. They derived a gauge field theoretical solution to 
the problem of a uniformly moving screw dislocation. However, we will show that their 
solution is not the correct one for a supersonic screw d is location. First of all, they u sed the 
constitutive relations given by IKadic fc Edelenl (11983T): lEdelen fc Lagoudad (|1988T). which 
are t oo simple. Even in Edelen's version (IKadic k Edelenl . Il983l : lEdelen fc Lagoudad . 
19881 ) of the gauge theory of dislocations there are at least three characteristic veloci- 
ties. For the anti-plane strain problem of a screw dislocation Edelen's theory possesses 
two characteristic velocities, namely the velocity of shear waves and a gauge theoreti- 
cal characteristic velocity. This makes it impossible without further simplifications and 
arguments to use a Lorentz transformation or at least it is not possib le to transform a 
movin g screw dislocation into a steady-state dislocation. Nevertheless, ISharma fc Zhang 
(120061 ) have constructed a solution of a screw dislocation using the Lorentz transformation 
with respect to a gauge theoretical characteristic velocity. In addition, they have consid- 
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ered the gauge theoretical characteristic velocity as an upper limit velocity. Finally, they 
solved the equations of motion of a screw dislocation in a steady state. The stress field of 
their solution does not possess the correct far field behavior because their solution is only 
given in terms of the characteristic gauge velocity and not in terms of the shear speed 
of sound . Neither no Mach cones appear in the supersonic regime. At best, the solution 
given by ISharma fc Zhangl (120061 ) is only valid in the subsonic regime if the characteristic 
gauge velocity is equa l to the shear speed of s ound. 



In the meantime, lLazar fc Anastassiadid (120081 . 120091 ) have presented an improved 
translational gauge theoretical formulation of dislocations. They have used the correct 
isotropic constitutive relations in the gauge theory o f disl ocations which are incomp lete 
and too simple in the books of iKadic fc Edelenl (119831 ) and lEdelen &: Lagoudasl (119881 ). In 
the present paper we derive the equations of motion of a screw dislo cation and we will solve 
them for a uniformly moving screw dislocation using the theory of lLazar &: Anastassiadis 
( 120081 ). We will show that a dislocation is a massive field in the gauge theoretical formu- 
lation. We will calculate subsonic and supersonic solutions of a uniformly moving screw 
dislocation. 



2 Gauge theory of dislocations 



In di s location dynamics we have the following state quantities of dislocation^ (jKosevich . 
19791 : lLandau & Lifschit j . Il989h 



Tijk ftikj ftij,k i 



(2.1) 



called the dislocation density tensor and the dislocation current tensor, respectively. These 
are the kinematical quantities of dislocations and they are given in terms of the incom- 
patible elastic distortion tensor and incompatible physical velocity of the material 
continuum V{. The dislocation density tensor describes the location and the shape of 
the dislocation core. The dislocation current tensor describes the movement of the 
dislocation core and it contains rate terms (%j and The dislocation density and 

the dislocation current tensors have to satisfy the translational Bianchi identities (square 
brackets indicate skewsymmetrization) 



ijk,l 



0. 



0. 



(2.2) 



The first equation means that dislocations do not have sources and the second one repre- 
sents that the circulation of the dislocation current is proportional to the time-derivative 
of the dislocation density. On the other hand, (12.21) are compatibility conditions ensuring 
that Tij k and 1^ can be given in terms of and Vi according to (12. ip . 

In the dynamical translational gauge theory of dislocations the Lagrangian is of the 
bilinear form (linear theory) 



C = T-W 



■PiVi 



2 Dijlij 



1 We use the usual notations: /3y t := dkfiij and /3, 



J%3 ■- 



1 

~ 2 



(2.3) 
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Here, the canonical conjugate quantities (response quantities) are defined by 



_oc _ oc_ _oc_ __ 9 J^L (oa\ 

<)>■; *«- o >• , • / " ' oi,/ Hijk - -or,,,: {ZA) 

where pi, a^, Iij, and H^k are the momentum vector, the force stress tensor, the disloca- 
tion momentum flux tensor, and the pseudomoment stress tensor, respectively. 

The Euler-Lagrange equations derived from the total Lagrangian C = C(vi, /3y, iy, T^) 
are given by 

+ <rf^ (2 - 5) 



We add to £ a so-called null Lagrangian, Cn = o'ijPij — Vi v h with the "background" 
stress of- and momentum p° as external source fields. Written in terms of the canonical 
conjugate quantities (I2.4p . Eqs. (I2.5P and (I2.6P take then the form 

Dyj + Pi = p^, (momentum balance of dislocations) (2.7) 
L>ij + Hij k k + aij = o°- (stress balance of dislocations). (2.8) 

Eqs. (12.71) and (12. 8p represent the dynamical equations for the balance of dislocations. 
Eq. (12 .7p is the momentum balance law of dislocations, where the physical momentum is 
the source of the dislocation momentum flux. Eq. (12.81) represents the stress balance of 
dislocations. Thus, the force stress and the time derivative of the dislocation momentum 
flux are the sources of the pseudomoment stress. The conservation of linear momentum 
appears as an integrability condition for the balance of dislocation equations. This can 
be seen by applying dt on (12. 7p and Oj on (12.81) and subtracting the last from the first one 

Pi ~ = 0, (force balance of elasticity), (2.9) 

where the time-derivative of the physical momentum is the source of the force stress. 

The linear constitutive relations for the momentum, the asymmetric force stress, the 
dislocation momentum flux tensor and the pseudomome nt stress of an isotropic and cen- 



trosymmetric medium are (ILazar fc Anastassiadisl . 120081 ) 



Pi = pvi, (2.10) 

Oy = \5ij(3 kk + n(pij + (3ji) + 7(/3y - pji), (2.11) 

Dij = di5ijl kk + d 2 {Iij + Iji) + d 3 (Iij - Iji), (2.12) 

Hij k = c\Tijk + ciiTjki + Tkij) + cs^SijTuk + SikTiji), (2-13) 

where p is the mass density and with 9 material constants p,, A, 7, C\, . . . , C3 and d\, . . . ,d 3 . 
Here p, X are the so-called Lame coefficients, 7 denotes the modulus of rotation, ci, . . . , C3 
are higher-order stiffness parameters and di, . . . , d$ are related to higher-order inertia due 
to dislocations. 
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The requirement of non-negativity of the energy (material stability) E = T + W > 
leads to the conditions of positive semi-defmiteness of the material constants. Thus, the 



material parameters have to fulfill the following conditions (ILazar fc Anastassiadid . 120081 ) 



d 2 >0, d 3 >0, 3d 1 + 2d 2 >0, (2.14) 

/i>0, 7>0, 3A + 2^>0, (2.15) 

ci-c 2 >0, ci + 2c 2 >0, ci-c 2 + 2c 3 >0. (2.16) 

Substituting the constitutive equations in the above system (12.71) and (12.81) . we obtain 

di($jj,i ~ v jd i) + (d 2 + d 3 )(/3 ijd - Vi dj ) + (d 2 - d 3 )( - v jdi ) + Pi = p°, (2.17) 
diSij(p kk - v kjk ) + (d 2 + d 3 )(Pij - Vij) + (d 2 - d 3 )((3. 



jl Vj,i 



+ Ci(PikJk — Pij,kh) + C 2 (f3ji^k — Pjk,ik + Pkj,ik ~ Pkijk) + C 3 [5ij(Plk,lk ~ Pll,kk) 

+ (Pkkji — Pkj,ki)] + a ij = a iji (2.18) 

which is a coupled system of partial differential equations for the field quantities v and 
P- 

3 Equations of motion of a screw dislocation 

We now proceed to derive the equations of motion for a moving screw dislocation. The 
symmetry of a screw dislocation leaves only the following non-vanishing components of 
the physical velocity vector and elastic distortion tensor: v z , f3 zx , (3 zy . The equations of 
motion of a moving screw dislocation read 





(d 2 + d 3 )($ zx> 


x + fizy,y 


- Av z ) + pv z 




(3.1) 


(d 2 + d 3 )((3 zx 


^z,x) + Cl(Pzy,xy 


~ Pzx,yyj 


+ (/i + i)(3 zx 




(3.2) 


(d 2 - d 3 )((3 zx 


^z,x) + C 2 (ftzx,yy 


~ Pzy,xy) 


+ l)f3zx 




(3.3) 


(d 2 + d 3 )((3 zy 


^z,y) + C\(y(5 ZXX y 


Pzy,xx , 


+ + j)(3 zy 


= in + 


(3.4) 


(d 2 - d 3 )0 zy 


^z,y) + £ 2 (ftzy,xx 


Pzx,xy ) 


+ ((l~ i)$zy 




(3.5) 



where A = d xx + d yy . In addition, the equilibrium condition is given by 

{(I + j)(Pzx,x + Pzy,y) = PV Z . (3.6) 

From the system of equations fl3.2l) - fl3.5p we obtain the two relations 

C\ c 2 



(i + l /i-7 
d 2 + d 3 _ d 2 - d 3 
/i + 7 /i-7 ' 



(3.7) 
(3.8) 
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We may introduce the following quantities 



f 2 - 01 

/i + 7 


(3.9) 


Ll ~ P ' 

2 G?2 + 
T l — • 

/i + 7 


(3.10) 


(3.11) 



Here l\ and L\ are the 'static' and 'dynamic' characteristic length scales and T\ is the 
characteristic time scale of the anti-plane strain problem. Moreover, l\ is related to 
dislocation stiffness and L\ is related to dislocation inertia^). The velocity of elastic shear 
waves is defined in terms of the 'dynamic' length scale Li and the time scale T\\ 



L\ /i + 7 



(3.12) 



Due to the presence of 7 the velocity of elastic shear waves is greater than in a theory 
with symmetric force stresses. Moreover , the y elocity of shear waves has a similar form in 
micropolar elasticity (see, e.g., Nowacki ( 19861 )) where the force stress is also asymmetric. 
In a similar way, we may introduce the following transversal gauge-theoretical velocity 
defined in terms of l\ and t\\ 



t? d 2 + d 3 



and therefore 



Li 



(3.13) 



(3.14) 



In the case 7 = we recover symmetric force stresses and from relations (13.71) and 
Q3.8fl : Ci = — c 2 and d 3 = 0. It can be seen that c\ = — c 2 gives the inequalities (12.161) 
with ci > and —c\ > 0. Thus, for nonnegative ci and c\ 7^ we cannot fulfill ( 12.161) . 
This is the price we have to pay for symmetric force stresses of a screw dislocation in the 
dislocation gauge theory. 



B ecause we deal with the phy sical state quantities (v, (3), no pseudo-Lorentz gauge (IKadic fc Edelenl . 



19831 ; lEdelen fc Lagoudad . Il988l ) is used and allowed during the simplification of the equa- 
tions of motion. Gauge conditions are only allowed for gauge potentials and not for 
physical state quantities. Of course, for anti-plane strain the equilibrium condition ( 13.61) 
together with ( 13.121) looks like a 'gauge' condition but it is not a gauge condition from 
the physical interpretation. 



2 We want to mention that similar length scales (£ s - static ch a racteristic length, £4 ~ dynamic char- 
acteri stic length) have been also obtained bv lGitman et al. (2005); lAskes fc Aifantisl (|2006l) ; lAskes et al 
(|20071 ) in a dynamic gradient elasticity. 
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Applying the equilibrium condition (13.61) . the equations of motion (I3.1l) -( l3~5l) can be 
written in the form 



t\v z - L'(Av z + v z 



-HI 



IA(3 ZX - r\ (l - ^jv ZjX + (3 ZX = (3° zx , 
lA(3 zy -r^{l-^)v z , y + (3 zy = f3° zy . 



(3.15) 
(3.16) 

(3.17) 



These are the equations of motion of an arbitrary moving screw dislocation in the frame- 
work of dislocation gauge theory. Some important questions arise if the lengths l\ and 
Li are independent or not and is there a physical reason to decouple the field equa- 
tions fl3TT5l) - fl3Tn) . 

Consider a screw dislocation moving in the x-direction. If we want to construct a 
solution which is consistent with the classical solution, we have to fulfill the condition: 



0, because only the classical dislocation density T° and dislocation 



current 1° are non-zero. If we use the field equations (13 . 1 5|) and (13.161) . we obtain: 



r 02 



±d tt -L\A + l 



LI 



-T 



To guarantee that (13. 18[) is fulfilled, we choose 



and 



d tt - t\A + 1 



(3.18) 



(3.19) 



Equation (13.19a ) ensures that I zx = and (13.19b ) with (I3.14p gives the relation ax = ct, 
which means that only one characteristic velocity ct survives for a screw dislocation in 
the gauge theory of dislocations. Thus, for a physically consistent solution we obtain the 
uncoupled Klein-Gordon equations 



J ZX1 



[l+£lD T ]v z -- 

[l+£lD T ](3 zy = /3° zy , 
with the following (1 + 2) -dimensional d'Alembert operator (wave operator) 



a. 



d tt - A. 



(3.20) 
(3.21) 
(3.22) 



(3.23) 



-T 



Moreover, the uncoupled system (13 .201) — (13 .221) with only one characteristic velocity pos- 
sesses a Lorentz symmetry. However, ct is not a limit velocity unlike the speed of light in 
spec ial relativ i ty. In field theories, Klein-Gordon equations describe massive fields (see, 
e.g., iRubakovl (120021 )). Thus, a dislocation is a massive gauge field. From the condition 
L\ — £i, we find for the inertia term of a screw dislocation 



Cl 



G?2 + d 3 = ~y = p£x 



(3.24) 



r 
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that it is given in terms of the characteristic length scale l\. Under these assumptions, 
the dynamical dislocation gauge theory of a screw dislocation possesses only one internal 
length scale l\. 



If we m ultiply Eqs. ( 13 . 20|) — (I3.22[) with Dt and using the 'classical' result (IGunther . 



19681 . 119691 ). we obtain 



[i + eln T ]n TVl 

[l + ^Dr]DrA 



zy,y> 



t: 





zxy,yi 



zxy,x 1 2 zy 



(3.25) 
(3.26) 

(3.27) 



as a set of fourth-order partial differential equations. As source terms only the classical 
dislocation density Tf! and dislocation current 7° ar e acting. Eqs. (I3.25I)-(I3.27I) have 
the two-dimensional form of Bo pp- Podolsky equations ( Bopp . 1940 ; Podolsky . 19421 ) (see 
also Iwanenko fc Sokolow in generalized electrodynamics, introduced by Bopp and 

Podolsky in order to avoid singularities in electrodynamics. 

In general, the velocity V of a screw dislocation might be subsonic or supersonic 
relative to the material space. A subsonic velocity lies in the range: < V < ct and for 
supersonic screw dislocations the velocity reads: ct < V. 



4 Uniformly moving screw dislocation 

We now study a screw dislocation moving with a uniform velocity V in the x-direction. 
Here V denotes the dislocation velocity in the material space. The dislocation moves 
relative to the material space, which plays the role of an 'aether'. If a screw dislocation is 
moving with the velocity V, then v z is the physical velocity of the material space, where 
the dislocation lives in, in order to transport the dislocation core to another position 
in the material space. Let x' be the coordinate in the direction of motion in the moving 
coordinate system. For the uniformly moving screw dislocation, then we use the transform 



with 



and we obtain the equations 



x' = x-Vt (4.1) 



d t = -Vd x , (4.2) 



and 



[1 - ^((1 - M*)d X , X> + d yy )] Vz = V Z , (4.3) 

[1 - - M*)d x>x , + Oyy ) J (3 ZX = (3° zx , (4.4) 

[1 _ - M 2 T )d x , x , + d yy )]p zy = P° zy , (4.5) 

[1 - £l((l - M 2 T )d x , x , + d yy )]T zxy = T Q zxy) (4.6) 

[1 - £l((l - M 2 T )d x , x> + d yy )]l zy = 1% (4.7) 
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with the Mach number of a moving screw dislocation relative to ct'- 

V 

M T = — . (4.8) 
ct 

Eqs. (I4.3I) - (I4.7I) are two-dimensional modified inhomogeneous Helmholtz equations (sub- 
sonic case) . In the supersonic case they change to one-dimensional inhomogeneous Klein- 
Gordon equations (see below). 



4.1 Subsonic case 



We now study a screw dislocation moving with a uniform ve locity V < ct in the x- 
direction (My < 1). The 'background' solution reads (see, e.g.. iGiintherl (119681 )) 



v° z = Vd y F°, 



6 C 

Hz 



zxy 



b5(x-Vt)8(y), 



-d F° 

u y r , 



fa = (3 2 d x F c 



-bV5(x-Vt)5(y), 



with 



F° = — -3-lnr T 



2tT# 



T 



and 



,-,.= v , {x _ vty + (3 2 y\ p T = 
In the subsonic region, the field equations (14. 3 1) - (14. 7\i are: 



1 - M 2 . 



[i-e 2 Md x , xl + d, 

[l-£ 2 1 (fad x , x , + d yy )](3 zx = {3° z: 
[l-t\{l3 2 T d xlx , + d yy )]p zy = t3\ 



ZX1 



zyi 



and 



[l-e^d^ + dyy)]^ 

[l-£l((3 2 d x , xl + d yy )]L 



zxy 



T. 



zxy 



(4.9) 



(4.10) 



(4.11) 



(4.12) 
(4.13) 
(4.14) 



(4.15) 
(4.16) 



Eqs. ( 14 . 1 2 [) — (14 . 1 6 j) can be easily solved with the inhomogeneous parts f!4.9[) using Fourier 
transform or other techniques. The solutions for the dislocation density and the disloca- 
tion flux of a screw dislocation are given by 



- zxy 



<-zy 



27r£ 2 /3 T 



b V 



(4.17) 
(4.18) 



Here, the symbol K n stands for the modified Bessel function of second kind (McDonald 
function) and of order n. In Fig. [1] the dislocation density is plotted for different speeds 
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(a) 



x 



Vt 



(b) 



x 



Vt 



Figure 1: Contourplot of the dislocation density T zxy of a screw dislocation moving with 
subsonic speed: (a) V = O.lcx, (b) V = 0.9ct- 



subsonic with respect to ct- The field of the dislocation density suffers a contraction 
when the value of its velocity approaches the velocity op- The curve of T zxy is a circle at 
V = and is a generalized ellipse at any other velocity V < op- The field is dilated in the 
directions orthogonal to the direction of motion and contracted along the line of motion. 
The solution for the elastic velocity is given by 



b (3 T Vy 
2tt rl 



r T 



T ' J 



(4.19) 



In Figs. [2] and [3] the physical velocity of a screw dislocations is plotted for subsonic 
velocities. It shows again a Fitzgerald contraction. It does not have a singularity. The 
physical velocity possesses extremum values depending on the dislocation velocity near 
the dislocation line. The solution of the elastic distortion reads 



l - 







b T y 

2tt r\ 
b p T (x - Vt) 
2^ 



h0 



"T 



T 



-K x 



(4.20) 
(4.21) 



All the elastic fields are nonsingular in the gauge theoretical framework. 

On t he other hand, the solutio n f)4.20p and (I4.2ip has the same form as the expression 
given by ISharma &: Zhang! (120061 ) if ct = clt and 7 = 0. 



4.2 Supersonic case 

We now consider the supersonic case: V > ct, (the supersonic case: Mr > 1)- Therefore, 
the term 0j, = — 7^ alters where 



7t 



1. 



(4.22) 



10 



, . -10 -5 5 10 . . -10 -5 5 10 

(a) x-Vt (b) x-Vt 

Figure 2: Contourplot of the elastic velocity v z of a screw dislocation moving with subsonic 
speed: (a) V = 0.1c T , (b) V = 0.9c T . 



Then the field equation for the dislocation density changes to 

[1 + £?(7t£W - d yy )]T zxy = T 2 %. (4.23) 

It is an inhomogeneous one-dimensional Klein-Gordon equation. The dislocation den- 
sity of a supersonic screw dislocat ion has the form of the corr esponding Green function. 



Therefore, it is given by (see, e.g.. Ilwanenko fc Sokolowl (119531 )) 



T zxy = ^— Jo ( ^ r 1 ^- ) H((Vt -x)- lT \y\). (4.24) 



Here, the symbol J n denotes the Bessel function of first kind and of order n and H is 
the Heaviside step function. It is non-zero just for Vt — x > Jt\v\- It builds a shear- 
wave Mach cone with the angle sin#r = fr- The dislocation density has a maximum of 
T zxy = b/[2 , jT^i\ on the Mach cone. Inside the Mach cone it oscillates with decreasing 
amplitude. For the dislocation current density we obtain I zy = —VT zxy with (14.241) . 
The visualization of the Mach cone of the dislocation density T zxy produced by the screw 
dislocation in the supersonic regime is plotted in Fig. |U 

The equations of the elastic fields, altering their character from elliptic to hyperbolic, 

are 

[1 + ^(tt^v - v z = v° z , (4.25) 
[1 + l\{lla x , x , - d yy )]f3 zx = f3° zx , (4.26) 
[1 + t\{l 2 T a x , x , - d yy )](3 zy = (3° zy , (4.27) 
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(a) 




(b) 



Figure 3: Three-dimensional plots of the physical velocity v z in units of bV /2tt£ 2 i3t'- 
(a) V = O.lcr, (b) V = 0.9c r . 



with the inhomogeneous parts (1A.4j) - (1A.6j) . The corresponding solutions read 



bV 



Ji 



y/(Vt - xY - j 2 y 2 



fizx 
fizy 



2h ^(yt-xf -ily 2 ~\ £i1t 

b y T ( V '{Vt - x) 2 - ~i\y 2 



2£i yftyi - x) 2 - 7^2 
b {Vt - x) 

2fi y/(Vt - x) 2 - -ily 2 



■h 
Ji 



a/(W-x) 2 - 7 ^ 2 



H((Vt-x)- lT \y\), (4.28) 
H((Vt-x)- lT \y\), (4.29) 
H{{Vt-x)- lT \y\). (4.30) 



It can be seen that (14.28p is non-zero just for Vt — x > 7t|?/|- Also it builds a shear-wave 
Mach cone with the angle sin 8t = y 1 - It is important to note that the classical singularity 
as Dirac delta function of supersonic dislocations on the Mach cone does not appear in 
our gauge theoretical result. The elastic distortions and the velocity have a maximum 
value on the Mach cone: v z = &Vy/[4Lf7r], (3 ZX = —by/[Al\^T]i v z and [3 ZX are also zero at 
y = 0. Inside the Mach cone they oscillate with decreasing amplitude. The visualization 
of the Mach cone of the physical velocity v z produced by the screw dislocation in the 
supersonic regime is plotted in Fig. [5j 



Let us mention that the solution of ISharma &: Zhangl (120061) does not possess th e cor- 
rect behavior in the supersonic region. The solution given by Sharma fc Zhangl (12006) does 
not show Mach cones, which are predicted in computer s i mula tions (jGumbsch fc Gao , 



1999 



Koizumi et al. 



tally (INosenko et al 



2007h . 



20021 : iLi fc Shi 12002 : iTsuzuki et all 120081 ) and found experimen- 



5 Conclusion 

In this paper, we have investigated a moving screw dislocation in the gauge field theory of 
dislocations. We derived the equations of motion of an arbitrary moving screw dislocation 
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x-Vt 

Figure 4: Contourplot of the dislocation density T zxy of a screw dislocation moving with 
supersonic speed V = V2ct- 




-20 -15 -10 -5 



x-Vt 

Figure 5: Contourplot of the physical velocity v z of a screw dislocation moving with 
supersonic speed V = V2c T . 
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in this field-theoretical framework. First, a coupled system of inhomogeneous Klein- 
Gordon equations is obtained in the dynamical case with two characteristic velocities. We 
have found one dynamical characteristic length scale Li, one static characteristic length 
scale l\ and one characteristic time scale T\. Later we have decoupled the field equations, 
following physical arguments, to construct a consistent solution. Due to these arguments 
we found L x = t x and = cy. The elastic fields v z , B zx and 8 zy and the fields of the 
dislocation core T zxy and I zy have the characteristic speed op- For a uniformly moving 
screw dislocation we have given analytical solutions for the subsonic and supersonic cases. 
For the supersonic case we found one Mach cone for the velocities ct- 
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A Supersonic screw dislocation in elasticity 

In the elasticity, a supersonic screw dislocation moves with the velocity: V > cj<. In 
symmetric elasticity the speed of shear waves reads: = p/p and in asymmetric elasticity 
it is given by = (p+j)/ p- In asymmetric elasticity the antisymmetric part of the stress 
tensor produces body couples. T he field equations for the elastic velocity and the elastic 



distortions read (jGuntherl . 119681 ) 



n T v° z = -bV d y 5(x - Vt)5(y), (A.l) 
n T f3° zx = bd y 5(x-Vt)6(y), (A.2) 
a T (3° zy = b 7 2 T dJ(x-Vt)5(y). (A.3) 

These are inhomogeneous wave equations describing massless fields. The supersonic so- 
lutions are 

(A.4) 
(A.5) 
(A.6) 

where 5 denotes the Dirac delta function. Thus, the classical solutions are shock waves 
produced by a supersonic screw dislocation. Both in front of and behind the shock front 
the material is undeformed and at rest. 



v° = 

u z 


bV y 


5(Vt- 


z-7tM), 




2 \y\ 


8° = 

yzx 


b y 
2M 


5(Vt- 




8° = 

r J zy 


b^T 
2 


8(Vt- 


v-7t\v\), 



References 

Askes, H. & Aifantis, E.C. 2006 Gradient elasticity theories in statics and dynamics - a 
unification of approaches. Int. J. Fract. 139, 297-304. 



14 



Askes, H., Bennett, T. & Aifantis, E.C. 2007 A new formulation and (^-implementation 
of dynamically consistent gradient elasticity. Int. J. Meth. Engng. 72, 111-126. 

Bopp, F. 1940 Die lineare Theorie des Elektrons. Ann. Phys. (Leipzig) 38, 345-384. 

Callias, C. & Markenscoff, X. 1980 The nonuniform motion of a supersonic dislocation. 
Q. Appl. Math. 38, 323-330. 

Edelen, D.G.B. & Lagoudas, D.C. 1988 Gauge theory and defects in solids, in: Mechanics 
and Physics of Discrete System. Vol. 1. North-Holland, Amsterdam. 

Edelen, D.G.B. 1996 A correct, globally defined solution of the screw dislocation problem 
in the gauge theory of defects. Int. J. Engng. Sci. 34, 81-86. 

Eshelby, J.D. 1956 Supersonic dislocations and dislocations in dispersive media. Proc. 
Phys. Soc. B 69, 1013-1019. 

Frank, F.C. 1949 On the equations of motion of crystal dislocations. Proc. Phys. Soc. (Lon- 
don) A 62, 131-134. 

Gao, J. 1999 An asymmetric theory of nonlocal elasticity - Part 1. Quasicontinuum theory. 
Int. J. Solids Struct. 36, 2947-2958. 

Gitman, I.M., Askes, H. & Aifantis, E.C. 2005 The representation volume size in static 
and dynamic micro-macro transitions. Int. J. Fract. 135, L3-L9. 

Gumbsch, R. & Gao, H. 1999 Dislocations faster than the speed of sound. Science 283, 
965-968. 

Giinther, H. 1968 Uberschallbewegung von Eigenspannungsquellen in der Kontinuums- 
theorie. Ann. Phys. (Leipzig) 21, 93-105. 

Giinther, H. 1969 Zur Theorie der Uberschall-Eigenspannungsquellen. Ann. Phys. 
(Leipzig) 24, 82-93. 

Giinther, H. 1988 On Lorentz symmetries in solids. Phys. stat. sol. (b) 149, 101-109. 

Hehl, F.W. & Kroner, E. 1965 Zum Materialgesetz eines elastische Mediums mit Mo- 
mentenspannungen. Z. Naturforschg. 20a, 336-350. 

Hirth, J.P & Lothe, J. 1968 Theory of Dislocations. McGraw-Hill, New York. 

Iwanenko, D. & Sokolow, A. 1953 Klassische Feldtheorie. Akademie-Verlag: Berlin. 

Kadic, A. & Edelen, D.G.B. 1983 A Gauge Theory of Dislocations and Disclinations. in: 
Lecture Notes in Physics. Vol. 174, Springer, Berlin. 

Koizumi, H., Kirchner, H.O.K. & Suzuki, T. 2002 Lattice wave emissions from a moving 
dislocation. Phys. Rev. B 65, 214104-1-9. 

Kosevich, A.N. 1979 Crystal Dislocations and the Theory of Elasticity, in: Dislocations 
in Solids Vol. 1. F.R.N. Nabarro, ed., North- Holland, pp. 33-165. 



15 



Kroner, E. 1968 Interrelations between Various Branches of Continuum Mechanics, in: 
Mechanics of Generalized Continua, IUTAM Symposium, E. Kroner, ed., Springer, 
Berlin, pp. 330-340. 

Kroner, E. 1976 Uber eine besondere Symmetrie zwischen Verschiebungen und Span- 
nungsfunktionen in der ebenen Elastostatik. Ing. Arch. 45, 217-221. 

Landau, L.D. & Lifschitz, E.M. 1989 Elastizitatstheorie, Band 7, Lehrbuch der Theoretis- 
chen Physik. Akademie-Verlag, Berlin. 

Lazar, M. & Anastassiadis, C. 2008 The gauge theory of dislocations: conservation and 
balance laws. Phil. Mag. 88, 1673-1699. 

Lazar, M. & Anastassiadis, C. 2009 The gauge theory of dislocations: static solutions of 
screw and edge dislocations. Phil. Mag. 89, 199-231. 

Li, Q. & Shi, S.-Q. 2002 Dislocations jumping over the sound barrier in tungsten. Appl. 
Phys. Lett. 80, 3069-3071. 

Markenscoff, X. & Huang, S. 2008 Analysis for a screw dislocation accelerating through 
the shear-wave speed barrier. J. Mech. Phys. Solids 65, 2225-2239. 

Nosenko, V., Zhdanov, S. & Morfill, G. 2007 Supersonic dislocations observed in a plasma 
crystal. Phys. Rev. Lett. 99, 025002-1-4. 

Nowacki, W. 1986 Theory of Asymmetric Elasticity. Pergamon Press, Oxford. 

Pan, Ke-lin & Chen, Zhi-da. 1994 A screw dislocation by nonlinear continuum mechanics. 
Applied Mathematics and Mechanics 15, 1093-1102. 

Podolsky, B. 1942 A generalized electrodynamics: part I - non-quantum. Phys. Rev. 62, 
68-71. 

Rubakov, V. 2002 Classical Theory of Gauge Fields. Princeton University Press, Prince- 
ton. 

Sharma, P. & Zhang, X. 2006 Gauge field theoretical solutions of a uniformily moving 
screw dislocation and admissibility of supersonic speeds. Phys. Lett. A 349, 170-176. 

Tsuzuki, H., Branicio, P.S. & Rino, J. P. 2008 Accelerating dislocations to transonic and 
supersonic speeds in anisotropic metals. Appl. Phys. Lett. 92, 191909-1-3. 

Weertman, J. 1967 Uniformly moving transonic and supersonic dislocations. J. Appl. 
Phys. 38, 5293-5301. 

Weertman, J. & Weertman, J.R. 1980 Moving dislocations. In: Dislocations in Solids, 
Vol. 3. F.R.N. Nabarro, ed., North-Holland, Amsterdam, pp. 1-59. 



16 



